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Several experimental studies have shown the presence of spatially inhomogeneous phase coexis- 
tence of superconducting and non superconducting domains in low dimensional organic supercon- 
ductors. The superconducting properties of these systems are found to be strongly dependent on 
the amount of disorder introduced in the sample regardless of its origin. The suppression of the 
superconducting transition temperature T c shows clear discrepancy with the result expected from 
the Abrikosov-Gor'kov law giving the behavior of T c with impurities. Based on the time dependent 
Ginzburg-Landau theory, we derive a model to account for the striking feature of T c in organic 
superconductors for different types of disorder by considering the segregated texture of the system. 
We show that the calculated T c quantitatively agrees with experiments. We also focus on the role 
of superconducting fluctuations on the upper critical fields H C 2 of layered superconductors showing 
slab structure where superconducting domains are sandwiched by non-superconducting regions. We 
found that H C 2 may be strongly enhanced by such fluctuations. 



INTRODUCTION 

Low dimensional organic conductors have attracted 
sustained interest regarding their rich phase diagram 
showing almost all known ground states with particu- 
larly a superconducting phase next to a magnetic or- 
dering [1-5], which is reminiscent of high-T c super- 
conductors. This proximity is believed to result in 
a phase segregation consisting of an inhomogeneous 
coexistence of electronic states such as the nanoscale 
charge stripes observed in copper oxide [6]. Meso- 
scopic scale inhomogeneities have been reported in the 
organic charge transfer salts (TMTSF) 2 X [7-9], so- 
called Bechgaard salts, and k(BEDT-TTF) 2 Y [3, 10] 
where TMTSF denotes the donor molecule tetramethyl- 
tetraselenafulvalene while BEDT-TTF stands for the 
bis(ethylencdithio)tetrathiafulvalene; X=C10 4 , PFg, 
Re0 4 ... and Y=Cu(NCS) 2 , Cu[N(CN) 2 ]Z (Z=Br and 
CI) are inorganic anions. These compounds are highly 
anisotropic: the former show marked quasi-one dimen- 
sional (qlD) electronic properties whereas the latter, ab- 
breviated as (ET) 2 , can be regarded as a layered (two- 
dimensional 2D) systems. The switching between the 
several low temperature electronic phases of these com- 
pounds can be controlled by applying hydrostatic pres- 
sure or by chemical substitution of cither the inorganic 
anion or the donor molecule. 

A close investigation of the boundary between the mag- 
netic phase and the superconducting one reveals a meso- 
scopic phase segregation where the system shows coex- 
istence of superconducting domains separated by mag- 
netic ones [5, 7-10]. Several experimental studies have 
reported that the interplay between these two phases 
is substantially dependent on the purity of the sample. 
Likewise, it has been found that this segregated struc- 



ture can be induced by varying the hydrostatic pressure 
around the critical value P c at which the magnetic phase 
collapses [7-11]. Disorder turns out to be a key pa- 
rameter for the stability of the superconducting state in 
these compounds. Despite the general consensus on the 
unconventional nature of the superconducting phase of 
different organic superconductors, the symmetry of the 
corresponding order parameter is still under debate [2]. 
Studying the effect of disorder on the superconducting 
phase may bring an answer to the puzzling question of 
the superconducting symmetry. In these molecular ma- 
terials, two types of disorder can be found: (i) Extrin- 
sic disorder resulting from the random imperfections ob- 
tained by introducing impurities or radiation damages, 
(ii) Intrinsic disorder related to the noncentrosymmet- 
ric character of the inorganic anions, such as ClOj in 
(TMTSF) 2 C104[12], or to the presence of internal de- 
grees of freedom such as partially disordered ethylene 
groups adopting several conformations in k(ET) 2 X [3]. 
These structural degrees of freedom, which are strongly 
dependent on the cooling history of the sample, have 
drastic effect on the superconducting state. A rapid cool- 
ing of the sample induces a strong depression of the su- 
perconducting transition temperature T c and may even 
destroy superconductivity as in (TMTSF) 2 C10 4 [2]. 

There have been several experimental studies deal- 
ing with the disorder in organic superconductors. In 
(TMTSF) 2 X, studies were focused particularly on the 
cooling rate induced disorder in (TMTSF) 2 C10 4 [13- 
15] and disorder generated by the anion substitution in 
(TMTSF) 2 (C10 4 )i- a; (Re04)x [16-19]. In k(ET) 2 X, dis- 
order was introduced by rapid cooling [20-23], chemical 
substitution of donor or anion [24-28] , X-ray and proton 
damage [28-30]. 

The superconducting transition temperature T c is 
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found to be reduced by increasing the amount of dis- 
order regardless of its origin [31, 32]. A theoretical anal- 
ysis [33] of the dependence of T c on disorder in different 
(ET) 2 X salts has shown that the suppression of T c obeys 
to the Abrikozov-Gor'kov (AG) law [34] which describes 
the suppression of T c in the presence of either magnetic 
impurities or non magnetic impurities in non-s-wave su- 
perconductors. However, for large disorder, clear depar- 
ture from the AG formula has been reported in several 
low dimensional organic superconductors [15, 22, 25, 30]. 
Alternative mechanisms have been suggested to account 
for this peculiar behavior [30, 35, 36]. A more detailed 
discussion is given below. 

The outcomes of the studies dealing with the effect 
of disorder on organic superconductors are consistent 
with the unconventional character of the superconduct- 
ing state since the latter is suppressed in presence of 
non-magnetic impurities [2]. However, these studies do 
not state clearly the symmetry of the superconducting 
order parameter. Experiments carried out on Bechgaard 
salts [8, 37, 38] show that the upper critical fields H c2 
along the most conducting axes a and b', H^ 2 and H^ 2 , 
are clearly times larger than the Pauli paramagnetic limit 
Hp above which U" 2 and H h c2 continue to increase with 
no sign to saturation. This behavior seems to rule out 
the singlet state and gives compelling indications of the 
presence of a spin triplet superconducting phase. How- 
ever, the question of the superconducting gap symmetry 
in Bechgaard salts is not settled yet. Despite the general 
tendency to accept the singlet character of the low field 
state, the symmetry of the high field phase continues to 
be a controversial issue [2-4, 37] . 

In this paper, we first complete the theoretical study 
presented in Ref. [36] dealing with the role of cooling in- 
duced disorder on the inhomogeneous superconducting 
phase of (TMTSF) 2 C10 4 where it is found that T c is 
strongly suppressed as the mixed character is enhanced. 
We then present a quantitative comparison with experi- 
ments which was lacking in Ref. [36] . In the present work, 
we give a detailed analysis of the experimental data in 
quasi-lD and 2D organic superconductors to infer the 
key parameters used in our model. We then derive the 
behavior of T c as a function of various types of disorder. 
To the best of our knowledge our results are the first to 
account quantitatively for the suppression of the super- 
conducting transition temperature by various sources of 
disorder in several organic superconductors. There is a 
general consensus, emerging from recent experiments, on 
the discrepancy between the experimental data obtained 
at large disorder rate and the results expected from the 
AG law [15, 22, 25, 30]. Our work will shed light on 
the origin of this discrepancy. We also briefly present in 
this paper, the effect of inhomogeneity on the upper crit- 
ical fields in the case of layered organic superconductors 
where superconducting domains are sandwiched between 
non-superconducting regions as it has been reported in 



(TMTSF) 2 PF 6 [8, 9]. In the next section, we summa- 
rize the experimental results related to the dependence 
of T c on disorder in several low dimensional organic su- 
perconductors. After a brief comment on the theoretical 
results presented in literature, we describe, in Section 
3, the method we used to extract, from the experimen- 
tal data, the theoretical parameters of the model we de- 
scribed in Ref. [36]. In section 4, results are discussed in 
connection with experiments. It is worth noting that, 
as done in Refs.[39, 40], we do not consider a particu- 
lar symmetry of the superconducting gap. However, the 
decrease of T c as a function of non magneticc disorder, re- 
ported in the present paper, supports the scenario of un- 
conventional superconducting order parameter in the dif- 
ferent low dimensional organic superconductors discussed 
in this work. This scenario seems to be the subject of a 
general consensus emerging from the huge experimental 
results dealing with these compounds [2, 3, 41, 42]. Nev- 
ertheless, the debate on the gap symmetry is not com- 
pletely settled. 

SOME EXPERIMENTAL FACTS 

When cooled sufficiently slowly, the C10 4 anions in 
(TMTSF) 2 C10 4 order for entropy reasons [5, 43] at a 
structural transition temperature Tao ~ 24 K. The 
nature of the electronic ground state, stabilized below 
Tao, is found to be substantially dependent on the 
cooling rate. In the relaxed state (R-state), obtained 
for a cooling rate less than O.lK/min, the ground state 
is superconducting below 1.2 K [1, 5]. However, the 
quenched state (Q-state) resulting of a cooling at a rate 
more than 50K/min, is an insulating spin density wave 
(SDW) [13, 14, 44, 45]. The most striking effect of the 
cooling rate is obtained for intermediate cooling rates at 
which the sample exhibits an inhomogeneous mixture 
of superconducting domains (in which CIO4 anions 
are ordered) and SDW regions (where CIO4 anions 
are disordered) [13-15, 17, 19, 46-52]. Meissner effect 
signal [48, 49] showed that, by increasing the cooling 
rate, the insulating regions get larger at the expense 
of the superconducting ones. These results have been 
corroborated by X-ray study [46] showing a variation 
of the fraction of ordered CIO4 domains as a function 
of the cooling rate. The fraction value was found to be 
identical to that of the superconducting volume deduced 
from the Meissner measurement [51]. The X-ray study 
has also provided the average size of these ordered 
regions as a function of the cooling rate. 

This Phase segregation has been also induced by 
X-ray irradiation [52] and chemical substitution[17, 
18] in (TMTSF) 2 C10 4 and by hydrostatic pressure in 
(TMTSF) 2 PF 4 6 [7-9]. 

It is worth to stress that the possibility of 
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a SDW/supcrconducting phase coexistence in 
(TMTSF) 2 PF 6 was first suggested by Greene et al. 
[53] who also considered a filamentary superconducting 
domain structure. 

The signature of such coexistence has been 
also reported in the organic superconductor 
(MDT-TS)(AuI 2 ) .44i, where MDT-TS denotes 
methylenediothio-tetraselenafulvalene, showing a 
temperature-pressure phase diagram reminiscent of 
those obtained in (TMTSF) 2 X and k(ET) 2 X salts [54]. 

Several experimental studies have revealed that the 
cooling rate has a drastic effect on the low temperature 
electronic properties of k(ET) 2 X (X=Cu[N(CN) 2 ]Br, 
Cu[N(CN) 2 ]Cl and Cu(NCS) 2 ) [3, 5, 20-23]. Different 
studies of the hydrogenated K(ET) 2 Cu[N(CN) 2 ]Br and 
the dcuterated analog K(d 8 -ET) 2 Cu[N(CN) 2 ]Br have 
given evidence of inhomogeneous superconducting state 
whose volume fraction is found to be reduced by a 
quenching the sample [20, 21, 24, 25, 55-58]. Such 
inhomogeneity may explain the anomalous behavior 
of the inplane conductivity in these materials [59]. 
It has been also found that the inhomogeneous su- 
perconductivity in k(ET) 2 X compounds can also be 
generated by chemical [24-26, 31, 32, 55, 60] or hydro- 
static pressure [10, 11, 61-63] and irradiation [28, 30, 64]. 

The outcome of the experimental studies on k(ET) 2 X 
and (TMTSF) 2 X materials is that the proximity of the 
superconducting state to a magnetic insulating phase is 
a key ingredient for the formation of an inhomogeneous 
phase. In the following, we present a theoretical approach 
to account for the dependence of the superconducting 
transition temperature on the disorder amount regard- 
less of its origin. We assume that the disorder does not 
introduce local magnetic impurities in the ordered do- 
mains which is consistent with experimental findings. 



THE MODEL 

In this section, we first give a brief review of the the- 
oretical approaches proposed in the literature to explain 
the behavior of T c as a function of disorder in organic 
conductors. We then present the outlines of the model 
we derived in Ref.[36] to qualitatively interpret the de- 
pendence of T c on the cooling rate in (TMTSF) 2 C10 4 . 
The key point lacking in Ref. [36] is the correspondence 
between the theoretical parameters and the experimen- 
tal disorder texture necessary to provide a quantitative 
comparison with the experimental results. In the remain- 
der of the section, we propose a method the extract the 
values of these parameters from the experimental data. 



Theoretical approaches: brief review 

Powell and McKcnzie [33] have analyzed, in the frame- 
work of the AG law, the experimental results dealing with 
the effect of different types of disorder on the supercon- 
ducting transition temperature T c of quasi-2D organic 
conductors. It is worthwhile reminding that the AG law 
describes the suppression of T c by magnetic impurities 
for s-wavc pairing. For non-s-wave superconductors, T c 
is suppressed by non-magnetic impurities according also 
to the AG formula given by [34]: 

ta S-*G)-*G + 5i5i?) (1) 

where T c q is the superconducting critical temperature 
in the pure limit and ip(x) is the digamma function, 
r is the quasi-particle lifetime due to scattering with 
magnetic impurities (r = tm) or non- magnetic impuri- 
ties (t = tjv). Powell and McKenzie [33] compared the 
experimental data obtained in /3(ET) 2 X and k(ET) 2 X 
with various types of disorder to the expectations of the 
AG formula. They found an excellent agreement using 
two free parameters to fit the theory. However, their 
study was not conclusive concerning the conventional 
nature of the superconducting order parameter, since 
the agreement with AG law holds for both magnetic 
and non-magnetic impurities. We emphasize that the 
data analyzed by Powell and McKenzie are in the clean 
limit corresponding to small amounts of disorder which 
acts as scattering points. It is worth noting that Joo et 
al. [18] have also reported that the suppression of T c 
with increasing impurity concentration x in the solid 
solution (TMTSF^ClO^i^ReO^ is consistent with 
AG formula. Moreover, based on EPR measurements, 
the authors [18] ruled out the magnetic nature of the 
scattering centers. The sensitivity of T c to the impurity 
concentration was, therefore, taken as a signature of the 
unconventional character of the superconducting gap in 
(TMTSF) 2 C10<i. Despite the good agreement between 
the AG law and the experimental results in the clean 
limit of organic conductors, a clear discrepancy has been 
reported in the dirty limit of quasi-2D salts [22, 25, 30]. 
To account for this discrepancy, Analytis et al. [30] 
suggested to use a generalized AG formula involving two 
order parameters with both s-wave and unconventional 
components. The presence of two gradients in the impu- 
rity induced suppression of T c was ascribed to this mixed 
superconducting ordering. The rapid decrease of T c for 
small disorder amount is attributed to the non-s-wave 
part whereas the slowing down of the T c suppression is 
assigned to the s-wave component. However, to the best 
of our knowledge, there is no experimental evidence of 
such mixed superconducting gap in organic conductors 
[33]. 
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Time Dependent Ginzburg-Landau approach 

Recently, we have proposed a model to describe the 
suppression of the superconducting transition tempera- 
ture with increasing cooling rate in (TMTSF) 2 C104 [36]. 
Basically, the model concerns a layered superconductor, 
where each layer, corresponding to the most conduct- 
ing plane, has a phase segregation structure with super- 
conducting domains embedded in a non-superconducting 
host. For simplicity, we assume that the supercon- 
ducting clusters are identical and form a square lat- 
tice. These superconducting regions are interacting via 
Josephson couplings parameterized by J\ and J 2 along 
the inplane directions a and b respectively. These do- 
mains are also coupled along the transverse direction c 
by the interplane Josephson coupling Jo. This struc- 
ture is consistent with the one proposed by Miiller et 
al. [63] from fluctuation spectroscopy measurements 
in k(ET) 2 Cu[N(CN) 2 ]C1. The authors have also ar- 
gued that the mixed state of superconducting and non- 
superconducting regions can be regarded as a network 
of Josephson coupled junctions connecting the supercon- 
ducting clusters. 

However, one should emphasize that in real samples 
the superconducting domains have not necessary pla- 
nar structure with identical shape and size. For exam- 
ple, X ray measurements in (TMTSF) 2 C10 4 [46] sug- 
gested ellipsoidal structure. One should also consider 
randomly distributed superconducting clusters instead of 
the square lattice assumed in our model. Actually, tak- 
ing into account the thickness of the superconducting do- 
main along the least conducting axis will not change the 
overall behavior of our results since superconductivity in 
the considered materials has essentially a two-dimension 
(2D) character. Moreover, the shape of the supercon- 
ducting domains is not a key factor in our model as far 
as the corresponding size is larger than the coherence 
length to avoid the phase fluctuations complexity. The 
basic parameter is the size of the non-superconducting 
junction on which depend the Josephson couplings. It 
should be interesting to study the effect of a random dis- 
tribution of the superconducting clusters which may be 
introduced by assuming a spatial weight dependence of 
the inplane Josephson coupling parameters. However, 
such correction is not expected to affect the outcomes 
of the present model. Indeed, the results depend on pa- 
rameters extracted from experimental data which can be 
regarded as mean values associated to the randomly dis- 
tributed superconducting islands in the real material. It 
turns out, that considering different sized superconduct- 
ing clusters with a random distribution may improve the 
quantitative agreement between our results and experi- 
ments but does not affect the general behavior. 

In (TMTSF) 2 C10 4 and k(ET) 2 X salts, the Josephson 
junction is basically antiferromagnetic. Josephson effects 



through magnetic junctions have been discussed by An- 
dersen et al. [65]. 

The anisotropic structure of the organic superconduc- 
tors is at the origin of strong superconducting fluctu- 
ations which are considerably enhanced as the dimen- 
sionality of the system is reduced [66] . Superconducting 
fluctuations, which can be brought out in many phys- 
ical properties, express the presence of superconduct- 
ing precursor effects in the normal state far from the 
superconducting transition temperature. These fluctua- 
tions smear out the sharp superconducting phase transi- 
tion anomaly observed in 3D isotropic superconductors. 
In anisotropic systems, this transition is broadened and 
spread out over a transient regime where superconduct- 
ing fluctuations take place [66]. Several experimental 
studies have given evidence for the presence of super- 
conducting fluctuations in organic conductors [3, 67, 68]. 

The critical superconducting fluctuations occur in the 
vicinity of T c over a thermal interval ST characterized by 
large fluctuations of the order parameter. The strength of 
the critical fluctuations is given by the so-called Ginzburg 
number [661: 
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where the anisotropic parameter 7 = 4jJ- is the ratio of the 
inplane coherence length over the out-of-plane one, A(0) 
is the zero temperature inplane penetration depth and 
</>o is the quantum flux. For conventional superconduc- 
tors Gi is about 10 -8 [5] whereas for hight-T c cuprates, 
d amounts to 1CT 1 [5]. For (ET) 2 X and (TMTSF) 2 X 
salts, Gi is of the order of 10~ 2 [3, 5]. Due to the super- 
conducting fluctuations, T c is reduced compared to the 
mean-field value Tq obtained in a bulky superconductor 
in the absence of fluctuations [66]. To can be estimated 
as T = T* X P + G t T^ xp where T^p is the experimental 
critical temperature obtained in the pure limit. 

The critical dynamics of the superconducting order pa- 
rameter is governed by the Time Dependent Ginzburg- 
Landau (TDGL) equation [40] 



r, 



dt 



dF 



"I" Cnij j i)j 



(3) 



where ipnij is the superconducting order parameter in 
the domain of the n th plane. Here i (j) labels 

the superconducting island along the a (b) direction. 

ksT is the relaxation rate of the order 
parameter, £0 is the inplane coherence length and m is 
the effective pair mass in the (ab) plane. The Langevin 
forces in Eq.3 govern the thermodynamic fluctua- 
tions of the superconducting order parameter. F is the 
free energy of the superconducting phase (F s ) compared 
to the normal state (F norm ) which is given by [36]: 
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The coefficients a and 6 are given by a = a^e and b = 
^QK 2 e^h 2 /2m 2 , where ao = h 2 /2m£ 2 and e is the reduced 
temperature e = ln(T/Tb), K = is the GL parameter, 

e = 2e is the pair electric charge. L\ (L 2 ) is the size of 
the superconducting domain along the a (b) axis. The 
Josephson couplings are 



Ji 



2m*d{ 



2m* ' 



and Jq = 



2m c s 2 
(5) 

Here d± (d 2 ) is the distance between two next neighbor- 
ing superconducting domains along the a (b) axis, while 
s is the interplane distance(Fig.l). m* and m c are the 
effective pair masses in, respectively, the superconduct- 
ing domain and along the transverse c axis. Hereafter, 
we consider m* = m since they correspond to the in- 
plane effective Cooper pair masses whereas m c will be 
expressed as a function of the anisotropy parameter 7 as: 
7 = W£c = y/m c /m. 




FIG. 1. Geometry of the inhomogeneous superconducting 
structure adopted in the model. The superconducting islands 
(denoted SC) are embedded in a non-superconducting ma- 
trix. We assume that the superconducting regions form an 
array for simplicity. The size of the superconducting domains 
is L\ (L2) in the a (b) direction. The thickness of the non- 
superconducting domain is denoted by di (tfe) along the a 
(b) axis. The interplane distance s has been exaggerated for 
clarity. The intraplane (Ji and J2) and the interplane (Jo) 
Josephson couplings are also represented. 

The Langevin forces are correlated through the 



Gaussian white- noise law [40] : 

(U(r,t)C^(f ',f)> = 2T 1 k B TS(r- f')S(t i')^ 

with r = (X + i(Li + di),Y + j(L 2 + d 2 ),ns) and f ' = 
X + i'(L! + di),Y + f(L 2 + d 2 ),n's). 

Assuming a Hartree approximation for the quartic 
term in Eq.4 leads to a linear problem where the 
b\ipnij\ 2 ' l pnij term in the derivative of the free energy is 
replaced by &(|^mj| 2 )^mj [40]. Given the Eqs. 3 and 4, 
the TDGL equation can be written as: 
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b(\lpni 3 \ 2 }lpmj + Jl {2lpnij ~ 4>ni+l j 



- 1pni-l j) 



J2 (2^nij 
Jq {2ll) n ij 
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(7) 



It is worth stressing that the TDGL theory we used 
was first proposed by Puica and Lang [40] to study the 
effect of the superconducting fluctuations on the inplane 
conductivity in a self consistent Hartree approximation 
and by Mishonov et <J.[69] within a Gaussian approx- 
imation. Moreover, the non-Ohmic effect of the high 
electric field on the transverse magnetoconductivity has 
been studied in the frame of TDGL by Puica and Lang 
[70]. The authors have also interpreted, within the same 
approach, the excess Hall conductivity in the layered 
superconductors for an arbitrarily inplane electric field 
[71]. The TDGL has proven to be a powerful tool to in- 
vestigate the role of superconducting fluctuations in the 
critical region to go beyond the usual Aslamazov-Larkin 
and Lawrence-Doniach approaches. 

The form of the free energy given by Eq.4 is based on 
the model proposed by Puica and Lang [40] to study the 
role of critical fluctuations on the conductivity of high-T c 
superconductors. The authors considered a layered struc- 
ture with only interplane Josephson coupling Jo connect- 
ing homogeneous superconducting layers. Following the 
approach of Puica and Lang [40], the critical tempera- 
ture T c of the superconducting state whose free energy 
is given by Eq.4 can be derived by solving the equation 



(8) 



To solve this equation, we use the Green function method 



[40] and define 
Cnij (r, t) by: 



the Fourier transforms of ip n ij and 
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i/> nij (x,y,t) = J ^2 J^^ q (k x ,k y ,t)e lk * x e e" 



Cqi^x, k y , t) 
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where k = (k x , k y ) and d = a + b(\ip n ij\ 2 }. 



of Eq.7. 



Equation 9 is obtained by taking the Fourier transform We define the Green function R q (k, t, k' x ,t') of Eq.9 by: 
I 



d h 2 k 2 

^dt + ~2m +ai+ 2Jl ^ ~ coa ^ Ll + 



R q (k, t, k' x , t') = 5{k x - k' x )S(t - t') 



(10) 



where 



The solution of Eq.9 is given by: 

i> q {k,t) = [ d£ [ dk' x R q {k,t,k' x J)t q {k' x ,k y ,t') 



tpq{k, t) can be derived by taking the Fourier transform 
of the Green function with respect to time [40] 



R q (k,oj,k' x ,t')= / dtR q (k,t,k' x ,t')e iu ^-^ 



h 2 k 2 

a\ = aH -+2J2II— cos(k v (L2+d2))+2Jo(l— cosfqs). Eq.10 can be written, then, as: 

2m 

I 



H 2 k 2 

-iusT^ 1 + + ai + 2Ji(l - cos(fc„(Li + di)) 



R q (k,t,k' x ,t') = 5(k x -k' x ) 



(11) 



I 

One can then obtain (\ip n ij\ 2 ) and solve Eq.8 for e = calculations give rise to the following equation to which 
to derive the transition temperature T c . Straightforward T c obeys: 

I 



ln ±c + zKBic9 dQ dE ^ / dw f w 2 + [ E + 2 J' (1-cosip) 
To 7TS J J J 

+ 2J[ (l - cos Vffcosflj + 2J' 2 (l - cos ^ d2 \/E 



= 



(12) 



where E c and w c arc cutoff parameters which depend on the energy scales of the system while g and the reduced 
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Joscphson couplings J[ [i = 0, 1, 2) are given by: 



9 



and .][ = —. 

a 



From Eq.12, one can derive the dependence of T c as a 
function of the J[ couplings (Fig.l of Rcf. [36]). However, 
a quantitative comparison with the experimental results 
requires the determination of the relationship between 
the Josephson couplings and the segregated structure 
induced by disorder. This point is discussed in the next 
section. 

Let us now turn to the effect of the phase segregation 
on the upper critical held of Bechgaard salts. We focus, 
for simplicity, on the behavior of the upper critical 
field along the most conducting axis (H C 2 \\ a) to avoid 

I 



complexity related to the field induced confinement 
occurring for H \\ b' . Moreover, we will tackle the 
case of (TMTFS) 2 PF 6 in the vicinity of the critical 
pressure P c at which the SDW phase vanishes [2]. In this 
region, the inhomogeneous structure is well established 
and the geometry of the domains has been deter- 
mined experimentally [9]. We give in the following, the 
outlines of the model. Details will be published elswhere. 



In the slab structure, reported in Refs. [8, 9], supercon- 
ducting domains of thickness L are sandwiched by SDW 
insulating regions along the most conducting axis a. We 
denote by d the size of these SDW domains. The free 
energy of the system in the presence of a magnetic field 
along the a direction H = (H, 0, 0) is given by: 



F = F s - F norm =Y, J dX J dy J dz a\i> n \ 2 + ^| (v - ^W) V>«| 2 + - V>»h 



ii 2 + dvg 4 



. (13) 



i 

where tp n is the superconducting order parameter in the Considering the gauge A = (0, 0, Hy), the TDGL equa- 
n th slab and J is the Josephson coupling constant be- tion becomes: 
tween two neighboring slabs: J = 2 m'd' 2 • 
I 



dlp n ti 2 



eo 



^n = <n{f,t) (14) 



Using the method introduced by Ullah and Dorsey [39] 
and by Puica and Lang [72] in the case of layered super- 
conductors where superconducting planes are coupled via 
Josephson parameters, one can deduce the upper critical 
field H^ 2 for a superconductor with a slab structure. De- 
tails of calculations will be given in a forthcoming paper. 
The results is discussed in the next. 



Lessons from experiments 

Disorder can be generated, as we have seen, by dif- 
ferent sources: cooling, anion or donor substitution and 
irradiation. Regardless of the origin of the introduced 
disorder, we have to establish the connection between the 
experimental parameters measuring its amount and the 
size of the non-superconducting domain on which depend 
the Josephson couplings J\ and J 2 given by Eq.5. For 
simplicity, we assume hereafter that the superconducting 
clusters are isotropic with a length L. This assumption 
does not affect the outcomes of the study as we shall show 



r 



in the following. The size of the non-superconducting do- 
main is denoted, henceforth, d. 



The case of (TMTSF) 2 CIO4. 

Pougct et al. [46] have carried out a high resolution 
structural study on (TMTSF) 2 C104 based on X-ray syn- 
chrotron radiation diffraction. The authors determined, 
for different cooling rates, the average size Li (i~ a, b, 
c) along the crystallographic directions and the fraction 
/ of the ordered CIO4 domains where superconductivity 
develops at low temperature. From these data, we have 
derived the typical size of the non-ordered domains as- 
suming that the ordered domains have a cubic shape with 
a length L separated by non-superconducting slabs of 
thickness d along the three directions. This 3D structure 
seems to be in contradiction with the assumption of a lay- 
ered superconductor considered in our model. Actually, 
the superconducting phase in organic conductors has a 
3D character but regarding the relatively small coherence 
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length £ c along the least conducting axis c, this phase is 
substantially dominated by the inplanc parameters. The 
coherence lengths in (TMTSF) 2 C104 along the crystal- 
lographic axes are typically of £ Q ~ 800A, ^ ~ 300A and 
£ c ~ 20A [1, 46]. 

£ c amounts to about the interplane distance c = ISA 
(denoted s in Eq.5) which supports the assumption of an 
interlayer Josephson coupling whatever the cooling rate. 

Considering the isotropic structure, the fraction of the 
CIO4 ordered domains can be written as: 



/ = 



L 



L + d 



(15) 



the values of d in order to have a quantitative compari- 
son between our theoretical results and experiments. In 
Table 1 we have listed, for different cooling rates, the 
experimental values of the superconducting volume frac- 
tion /, the lengths of the ordered domains along the crys- 
tallographic axes denoted L a , Lt and L c and the aver- 
age length L = (_L L;,_L C ) 3 of a superconducting domain 
within the isotropic assumption. We also give the typical 
size d of the non-superconducting domain and the renor- 

malized Josephson coupling , where J' = — = ( *r 



(Eqs.4,5). These values are calculated for an average in- 
plane coherence length £ = VCo^b = 500A. 



Given the experimental data for / and L, one can deduce 



Cooling rate (K/min) 


f 


L a (A) 


L b (A) 


L C (A) 


L = (L a L b L c )i (A) 


d(A) 


4- xio 2 

En 


J' 


0.5 


0.95 


700 


1050 


1150 


950 


16 


3.2 


976.5 


1 


0.9 


550 


900 


900 


750 


27 


5.4 


342.9 


3 


0.7 


400 


650 


550 


500 


63 


12.6 


62.9 


5 


0.48 


300 


550 


450 


420 


116 


23.2 


18.6 



TABLE 1: Superconducting fraction /, typical sizes Li and the corresponding average L value of a superconducting domain 
in (TMTSF)2C104 at different cooling rates. The interdomain distance d and the renormalized Josephson coupling J' are also 

indicated (data of Ref.[46]). 



In our model, the average size L of the superconduct- 
ing domain is assumed to be larger than the coherence 
length £0 to avoid any complexity related to phase fluctu- 
ations of the superconducting order parameter. Accord- 
ing to the data of table 1, this assumption breaks down 
for a cooling rate greater than 3 K/min for the average 
inplanc coherence length £0 ~ 500A. Actually, for such 
rates, the Josephson coupling is no more efficient regard- 
ing the relatively large size of the non-superconducting 
junction. It then reduces to a tunneling effect. At con- 
stant Josephson coupling, the superconducting transition 
temperature tends to saturation as far as the size L of the 
superconducting domain is larger than the inplane coher- 
ence length £o- Indeed, numerical results [73] have shown 
that, in this case, the superconducting transition tem- 
perature is practically independent of L. However, for 



rapid cooling rate, experimental data [46] show that the 
size d of the Josephson junction is only slightly changed 
whereas the superconducting average domain size L is 
strongly reduced giving rise to the collapse of the super- 
conducting transition temperature. In the present work, 
we will focus on the slow and intermediate cooling rates 
for which T c exhibit a clear departure from the AG law. 
We will then assume that L is roughly unchanged by 
cooling and remains larger than £o- Despite this simpli- 
fying assumption, we will show in the next section, that 
the obtained results are quantitatively in agreement with 
the experimental ones. 

The case of k(ET)-zX salts 
In table 2, we list the different symbols used in the text. 



The effect of the disorder induced by the cooling rate 
on k(ET) 2 X (X=Cu(NCS) 2 and Cu[[N(CN) 2 ]Br) has 
been investigated by Yoncyama et al. [22] by measur- 



ing the inplane penetration depth. The authors con- 
cluded that the local clean approximation is in good 
agreement with experimental results, indicating that dis- 
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Symbol 


Definition 


L 


Size of the superconducting domain 


d 


Size of the non-superconducting domain 


f 


Superconducting volume fraction 


h 


in-plane mean free path a 


h(°) 


in-plane mean free path in the cleanest limit 
(smallest disorder rate) 


Po 


residual resistivity at which T c deviates from AG law 


I* 

'II 


in-plane mean free path corresponding to pj$ [30] 


d(0) 


Size of the non-superconducting domain in the cleanest limit 


8 


Size increase of the non-superconducting domain by increasing disorder rate 

d = d(0) + S 



TABLE 2 : Definition of the different symbols used in the text. 
a I II is taken equal to the size L of the superconducting domain assuming that the latter is free of impurities. This assumption 
is based on the idea that we only consider the effect of inhomogeneous superconductivity induced by disorder and disregard 

any possible local impurities inside the superconducting island. 



order induced by fast cooling leads to a short inplane 
mean free path {ii. The authors have also obtained, from 
the data of Stalcup et al. [21], the behavior of the su- 
perconducting transition temperature T c with the Dingle 
temperature Trj which is an indicator of the sample pu- 
rity regarding its dependence on the mean free path Z||: 
T D = 2^r, where Vf is the Fermi velocity. It is worth 
stressing that the inplane mean free path measured from 
the intralayer component of the conductivity may be dif- 
ferent from those deduced from cyclotron resonance ex- 
periments or de Haas-van Alphen effect as discussed by 
Singleton in Ref. [74]. It has also been found that the 
measured Dingle temperature contains the effect of spa- 
tial inhomogeneities [74]. 

The domains, where the ethylene endgroups arc disor- 
dered by quenching, act as scattering centers according 
to the results of Yoneyama et al. [22]. The faster the 
cooling, the shorter the mean free path. By quench- 
ing, T c is reduced whereas the penetration depth in- 
creases, in particular for the Cu[[N(CN) 2 ]Br salt. One 
can then describe the system as a bulk superconductor 
where the number of defect centers increases with cool- 
ing rate. These defects are actually non-superconducting 
domains of a size d where the ethylene groups are ex- 
pected to be disordered. Actually, the idea attributing 
the disorder in k-(ET) 2 salts to the quenching ethylene 
endgroups disorder has been widely supported [3]. This 
idea is not, however, corroborated by the structural inves- 
tigation reported by Wolter et al. [75] who have argued 



that the configurationally disorder of ethylene endgroups 
is not the unique factor. The detailed structural study 
of Wolter et al. [75] is in agreement with the results of 
Strack et al. [76] showing that the contribution of the 
frozen-in disorder of the ethylene group is not substan- 
tial. 

It is worth noting that the defects induced by the 
configurational disorder of the ethylene endgroups 
cannot be considered as local impurities since the 
network of ff-bondings of the anions with the disordered 
terminal ethylene groups induce a displacive disorder 
involving several ET molecules [12]. From this picture, 
the mean free path l\\ can be taken as the size of the 
superconducting domain in our isotropic model. Let 
us denote by iii(0) the mean free path associated to 
the cleanest limit and d(0) the corresponding distance 
between neighboring superconducting domains. At 
a given cooling rate, i|| writes as Z|| = i||(0) — S and 
as a consequence d = d(0) + S due to the isotropic 
structure of the superconducting cluster we assumed 
in our model. The experimental values of l\\ and Z||(0) 
can be deduced from the Dingle temperature (Fig. 8 of 
Ref. [22]). By adjusting l\\(0) to obtain numerically the 
critical temperature T c for the slowest cooling rate, 
one can simply extract the experimental values of d 
for different cooling rates and derive the theoretical 
dependence of T c as a function of the cooling rate or the 
Dingle temperature. 

In table 3 we give the main parameters inferred from 
the experimental data of Ref. [22]. 



Yoneyama et al. [25] have studied the effect of deuter- 
ated ET molecule on the superconducting state of /t-[(hg- 



ET) 1 _ a; (d 8 -ET) a; ] 2 Cu[N(CN) 2 ]Br where h 8 -ET and d 8 - 
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Dingle Temperature Td (K) 


Mean free path l\\(A) 


d(A) 


2.04 


372.5 


7.5 


2.25 


338 


42 


2.5 


304 


76 



TABLE 3 : Some values of the inplane mean free path lu and the size d of the non-superconducting domains deduced from 
the experimental values of the Dingle temperature which have been obtained by Yoneyama et al. [22] from the data of 
Stalcup et al. [21]. We have assumed that the cleanest limit corresponds to the highest transition temperature which is 
ascribed, according to the data of Ref. [22], to Td — 2 K. The latter gives i||(0) ~ 380 A. 



ET correspond, respectively, to the fully hydrogenated 
and fully deuterated molecules. Moreover, Sasaki et al. 
[24] have obtained, based on scanning micro region in- 
frared spectroscopy, the superconducting volume frac- 
tion / in this salt as a function of the substitution ra- 
tio x. We can, then, deduce, for each concentration x, 
the size d of the non-superconducting domain as (Eq.15): 

f = ^j^j 3 — 1. According to our calculations [36, 73], 
the superconducting transition temperature T c is, prac- 
tically, unchanged by varying the size L of the supercon- 
ducting cluster. The latter should be, however, greater 
than the inplane coherence length £o to avoid, as we have 
previously discussed, any problem related to the phase 
coherence of the order parameter as it is the case in gran- 
ular superconductors [77]. 



X 


superconducting 


i xlO 2 


d8-ET concentration 


volume fraction / 




0.5 


1 





0.6 


0.9 


3.6 


0.7 


0.82 


6.8 


0.8 


0.68 


13.6 


0.9 


0.54 


23 


1 


0.44 


31.5 



TABLE 4 : Values of the superconducting volume fraction / and the ratio extracted from the experimental values of 
Ref. [24]. Here | = U\ 5 - 1. We have taken, as in Ref. [25], the cleanest limit at x = 0.5. 



We take L ~ 100A whatever the x value since for the 
considered k(ET) 2 X salts £ ~ 30A - 70A. Given the 
experimental values of / for different concentration x, 
one can, then, deduce d and establish the correspondence 
between the Josephson couplings of our model and the 
substitution ratio x to derive the theoretical behavior of 
T c as a function of x. It is worth stressing that we will 
consider the region of inhomogeneous superconductivity 
corresponding to x in the range 0.5 — 1 in the phase 
diagram of Yoneyama et al. [25]. The cleanest limit, 
corresponding to a homogeneous superconducting state, 
is ascribed to x — 0.5 (Fig. 5 of Ref. [25]). 

In table 4 we give the main parameters deduced from 
the experimental data of Sasaki et al. [24] . 



The T c dependence on the irradiation damage has 
been studied by Analytis et at [30] in k(ET) 2 Cu(NCS) 2 . 
The defect density was probed by the residual resistivity 
po which is proportional to the scattering rate r. We 
denote by Pq and IT, the critical value of po and the corre- 
sponding mean free path above which T c decs not follow 
a linear behavior [30] . Analytis et al. [30] have found that 

for pX, jt ~ 0.2 where £q ~ 70 A is the inplane coherence 

1 

length. This leads to the estimation IT, ~ 350A For a 



given irradiation amount, we can deduce the mean free 
path as 2 1| = since po is proportional to the inverse 
of the inplane mean free path lu. It is worth noting that 
the determination of the the inplane scattering rate or 
the inplane mean free path from the interplane residual 
resistivity has been used by Joo et al. [18, 19] in the case 
of the quasi-one dimensional organic superconductor 
(TMTSF) 2 C10 4 to fit, using the Abrikosov-Gorkov 
formula, their experimental data dealing with the de- 
pendence of the superconducting transition temperature 
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with the interplane residual resistivity. The result of Joo 
et al. was taken as an evidence of the unconventional 
character of the superconducting order parameter in 
(TMTSF) 2 C10 4 . Moreover, the extraction of the inplane 
scattering rate and the inplane mean free path from 
interlayer resistance to fit Abrikosov-Gorkov law has 
been also used by B. Powell and R. Mackenzie [33] 
in the case of several quasi-two dimensional organic 
superconductors . 



path in the clean limit corresponding to the highest T c 
obtained for p Q ~ 0.33 flcm. This gives Z||(0) ~ 6 IT,. 
The distance d between neighboring superconducting re- 
gions reads, then: d = d(Q) + 5, where S = l\\(0) — l\\ ~ 

d(0) is taken from the fit of the T c values in 



IT, 



6 



Po 
Po 



the cleanest limit. 



Following the method we introduced above, l\\ can be 
written as In = fii(0) — S, where iii(0) is the mean free 



Table 5 shows the main parameters inferred from the 
experimental data of Analytis et al. [30] . 



J 



po (flcm) 


Mean free path Z|| = ZJ ^ (in the unit of ZjJ) 


5 — 2|| (0) — Z|| (in the unit of Z|j) 


0.33 


6 





0.66 


3 


3 


1.5 


1.3 


4.7 


2 


1. 


5. 


2.6 


0.7 


5.3 


6 


0.3 


5.7 



TABLE 5 : Some values of the mean free path Zy and the non-superconducting size increase S in units of It corresponding to 
the value po of the residual resistivity at which the superconducting transition temperature T c deviates from the AG law. 
According to the data of Analytis et al. [30], po ~ 2S7cm. Taking a coherence length of £o ~ 70 A gives Z| ~ 350A. At the 
cleanest limit, corresponding to the highest T c , we obtain Zy(0) ~ 6Z|. 



RESULTS AND DISCUSSION 



In figure 2 we plot the superconducting transition tem- 
perature T c obtained from Eq.12 as a function of the 
cooling rate for (TMTSF) 2 C10 4 salt. The results are 
in good agreement with the experimental data of Mat- 
sunaga et al. [51] and Joo et al. [15]. As shown in figure 2, 
T c falls with increasing cooling rate but presents a sig- 
nificant slope change at a cooling rate of about 5 K/min 
above which the effect of the cooling rate on T c is re- 
duced. For slow cooling (less than 5 K/min), the tiny 
non-superconducting domains act as local defects which 
are known, according to the AG law, to have drastic 
effect on superconductivity. Since EPR measurements 
[14, 15] revealed the non-magnetic nature of the cooling 
induced disorder centers in (TMTSF) 2 C104, the sensitiv- 
ity of T c to the cooling rate is a signature of the non-s- 
wave character of the superconducting order parameter. 
By quenching, the non-superconducting regions get wider 
and can no more be considered as point defects. The de- 
crease of T c is slowed down and dees no more follow the 
AG formula. 



The tendency to saturation in this regime reflects 
the insensitivity of the Josephson coupling due to the 
increasing size of the Josephson junction between the 
superconducting domains. Actually, this saturating 
regime is expected to end as soon as the size of the 
superconducting domain becomes smaller than the 
inplane coherence length leading to the disappearance 
of the superconducting state. It is worth noting that, 
despite the simplifying assumption considered in our 
model, the obtained results show good agreement with 
experimental data. This indicates that actually, in the 
range of validity of the calculations, the cooling rate 
acts principally on the size of the non-superconducting 
domains. It turns out that the ratio -j^, on which depend 

the reduced Josephson coupling constants J' = in 
Eq.12, is the key parameter governing the dependence 
of the superconducting transition on the cooling rate. 

In figure 3, we depicted the dependence of T c on the 
Dingle temperature corresponding to different cooling 
rates in K-(ET) 2 Cu[N(CN) 2 ]Br. The good agreement be- 
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FIG. 2. Superconducting transition temperature as a function 
of the cooling rate in (TMTSF)2C104. Squares correspond to 
our results obtained from Eq.12 whereas the triangles and 
filled circles are, respectively, the data of Refs.[15, 51]. Lines 
are guide to eye. 



tween our results and the experimental data [22] supports 
the method used to connect the microscopic parameters 
of our model with the experimental amount of disor- 
der. It also reflects that the key ingredient is the size 
d of the non-superconducting domains which governs the 
strength of the Josephson coupling. 
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11.4 
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FIG. 3. Superconducting transition temperature as a function 
of Dingle temperature for K-(ET) 2 Cu[N(CN) 2 ]Br. Squares 
correspond to the results obtained from Eq.12 whereas the 
filled circles are the data of Ref. [22]. T5 denotes the value 
of the Dingle temperature at which the experimental data 
deviate from the AG law. 



In figure 3 we also give the fit of the experimental 
data with AG formula. At low disorder amount (small 
To), the experimental results are well described by 
the AG law, indicating that the non-superconducting 
domains are tiny enough and can be considered as 
defect points. In this regime, our results are also in 
agreement with the AG formula which means that one 
recovers the AG law as a limiting case corresponding to 
small disorder amount. However, for a large disorder 
amount, there is a clear discrepancy between the AG 
law and the experimental data since, in this case, the 
defect point picture dees not hold anymore regarding 
the increasing size of the non-superconducting domains. 
The critical Dingle temperature , marked by an arrow 
in Fig. 2, corresponds to a ratio |- ~ 0.2. For d > d* , 
the non-superconducting domains are large enough and 
the AG formula is no more justified. From the value of 
Tp, one can also deduce the reduced critical value of the 
mean free path If, : p ~ 0.17. 



Figure 4 (a) and (b) show the behavior of T c 
with the deuteration rate x in K-[(h8-ET)i_ a: (d8- 
ET%]2Cu[N(CN)2]Br for the relaxed and quenched sam- 
ples respectively. Our results are consistent with the data 
of Yoneyama et al. [25] . The calculated T c seems to have 
a decrease rate somewhat greater than the experimen- 
tal values, in particular for the quenched sample. This 
may be due to the simplified assumption concerning the 
isotropic character of the superconducting domains which 
we used to derive the length of the non-superconducting 
regions from the experimental data. With increasing x, 
the non-superconducting islands develop resulting in a 
reduced T c . However, the effect of deuteration induced 
disorder appears to be rather negligible. The suppression 
of T c was ascribed by Yoneyama et al. [25] to the pres- 
ence of magnetic impurities. The non-superconducting 
domains should be AF since the ground state of the ds- 
ET salt (x = 1) is AF with minor superconducting re- 
gions. 

The suppression of T c with irradiation in 
k(ET) 2 Cu(NCS) 2 is shown in Fig. 5 where the residual 
resistivity p measures the irradiation damage. A fit 
by the AG formula is indicated by the dashed line. 
Our results describe correctly the experimental data 
of Analytis et al. [30] and those of Sasaki et al. [28], 
which corroborates the idea that irradiation generate 
non-superconducting grains in the bulk of the sample 
and not only at the surface. Moreover, the departure 
from the AG law can be understood within the present 
model reduction of the Josephson tunneling due 

to increasing size of the non-superconducting regions. 
There is no need to assume a multicomponent supercon- 
ducting order parameter as proposed by Analytis et al. 
[30], especially since there is no experimental evidence 
for such scenario. 
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FIG. 4. Superconducting transition temperature as a func- 
tion of the deuteration concentration x for K-[(hs-ET)i_x(d8- 
ET) I ] 2 Cu[N(CN) 2 ]Br in the relaxed (a) and quenched (b) 
samples. Squares correspond to the results obtained from 
Eq.12 and the filled circles are the data of Ref.[25]. 



At the critical value p* at which the data show a 

deviation from the AG law, Analytis et al. [30] have 

found that - 0.2. This is practically the same value 
ii 

we obtain in the case of K-(ET)2Cu[N(CN)2]Br where 
~ 0.17. It seems that there is a critical ratio which 

ii 

limits the range of validity of the AG law in k-(ET)2X 
salts. More experimental data are needed to generalize 
this finding. 



The outcome of these results is that our model, 
despite its simplicity, provides a coherent interprcta- 
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FIG. 5. Superconducting transition temperature as a func- 
tion of the residual resistivity in k(ET)2Cu(NCS)2- Squares 
correspond to the results obtained from Eq.12 and the filled 
circles are the data of Ref. [301. 



tion of the role of disorder in layered inhomogeneous 
superconductors. The technique used to extract the 
leading microscopic parameter of the model related to 
the non-superconducting domain length, has proved to 
be reasonable regarding the good agreement between 
the calculations and the experimental findings. 

Figure 6 shows the behavior of the upper critical field 
along the a axis for different rates of disorder assuming 
the slab structure found in the case of (TMTSF) 2 PF 6 . 
The parameter d is the size of the non-superconducting 
slabs. In the case of large disorder amount, correspond- 
ing to a large d, H C 2 saturates at a value around the 
Pauli limit (~ 2.5 T). However, by reducing the disor- 
der concentration, which turns out to reduce d, H C 2 is 
greatly enhanced with a non saturating behavior as ob- 
served experimentally [8]. We ascribe this feature to the 
superconducting fluctuations which are enhanced by re- 
ducing d giving rise to a robust superconducting state. 
Our results suggest that the superconducting fluctua- 
tions are responsible of the non-saturating behavior of 
H c2 in (TMTSF)2PF 6 salt. This finding may shed light 
on the origin of the discrepancy between the H C 2 val- 
ues obtained by thermodynamic and transport measure- 
ments in (TMTSF)2C10 4 [37, 78, 79]. The former reveal 
a saturating upper critical fields whereas the latter show 
a non saturating behavior. Transport measurements are 
sensitive to superconducting fluctuations. According to 
our results, the divergent character of the upper critical 
fields reported by transport measurements is a signature 
of superconducting fluctuations. Regarding the non sat- 
urating behavior of the upper critical field, the high field 
phase has been ascribed to a triplet state whereas the 
low field phase is assumed to be a singlet phase [4]. 
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It is worth to note, that, as done by Ullah and Dorsey 
[39] and Puica and Lang [72], we do not consider in our 
calculations a particular symmetry of the superconduct- 
ing state. The model deals with unconventional super- 
conductors regarding the dependence of the critical tem- 
perature on disorder. However, according to our result, 
the high field phase could be a singlet state and the di- 
vergent character is nothing but the signature of the su- 
perconducting fluctuations. 




FIG. 6. Upper critical fields in layered superconductor with 
a slab structure. Calculations are done for (TMTSF)2PF6 
data, d is the thickness of the non superconducting domains 
and £o is the coherence length (fo ~ 500A) 



CONCLUDING REMARKS 

We have studied the effect of disorder, regardless of its 
origin, on the behavior of the superconducting critical 
temperature T c in the Bechgaard and k-(ET)2X salts. 
We have considered an inhomogeneous superconduct- 
ing phase where superconducting domains are embedded 
into a non-superconducting matrix. For small disorder 
amount the latter can be considered as defect point and 
the suppression of T c is found to be well described by the 
AG law. However, in the regime of large disorder rate, 
the non-superconducting regions get wider and the re- 
sults expected from AG formula are in clear discrepancy 
with experimental data. We have then proposed a model 
to describe the suppression of T c taking into account the 
inhomogeneous structure of the disordered material. We 



have shown that this texture gives rise to Josephson tun- 
neling which depends on the key parameter ^ where d 
is the size of the non-superconducting junction and £o is 
the inplane coherence length. We have also found that, 
below a critical value rj* of the ratio r] = j 2 -, where Z|| is 
the inplane mean free path, the AG law holds. The non- 
superconducting regions act as impurities in this case, rj* 
is found to be of the order of 0.2 in k (ET^X salts. 

However, above rf , the AG formula is no more reliable 
and one should consider sizeable non-superconducting 
domains through which Cooper pairs tunnel between su- 
perconducting islands. We have proposed a rather sim- 
ple method to extract the key microscopic parameters 
from the experimental data to establish the relationship 
between our calculations and the experimental data. A 
good agreement has been found for various sources of 
disorder despite the simplifying assumptions. The model 
could be improved by including the effect of a random 
distribution of the superconducting domains as observed 
recently in the high-T c superconductors La2Cu04 + y 
[80] where a structural ordering of the oxygen intersti- 
tials has been reported. This order, which is highly in- 
homogeneous, is characterized by a fractal distribution 
which seems to enhance the superconducting transition 
temperature [80]. Such effect should be also relevant in 
(TMTSF) 2 C104 where previous study [46] has shown the 
presence of a large distribution of the ordered CIO4 do- 
main size in the relaxed sample. It has been found that 
the width of this distribution shrinks as the cooling rate 
increases [46]. We have also studied the role of super- 
conducting fluctuations in layered superconductors with 
slab structure where transport measurements revealed 
non saturating upper critical field. We have ascribed this 
feature to the enhanced superconducting fluctuations to 
which transport probes are sensitive. Our results may ex- 
plain the absence of any divergent behavior in the upper 
critical fields reported by thermodynamic measurements 
in organic superconductors. 
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